Effective continuum threshold for vacuum-to-bound-state correlators 
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We formulate a new algorithm for obtaining the effective continuum threshold in vacuum-to- 
bound-state correlators — the basic objects for the calculation of hadron form factors in the method 
of light-cone sum rules in QCD. The effective continuum threshold is one of the key ingredients 
of the method which to a large extent determines the value of the form factor extracted from 
the relevant correlator. In a quantum-mechanical potential model, where the exact form factor is 
known, the application of our algorithm is shown to lead to a dramatic increase of the accuracy 
of the extracted form factor compared to the standard procedures adopted in the method of sum 
rules in QCD. Moreover, the application of our algorithm considerably enlarges the range of the 
momentum transfer where the form factor may be extracted from the correlator. 
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The extraction of the ground-state parameters from a correlator of quark currents is a cumbersome procedure: even 
if several terms of the expansion for the correlator (i.e., the OPE in the case of vacuum-to- vacuum correlators and 
the twist expansion in the case of vacuum-to-hadron correlators) are known precisely, the numerical procedures of the 
method of sum rules [11,01 cannot determine the true exact value of the bound-state parameter. Instead, the method 
should provide the band of values such that the true hadron parameter has a flat probability distribution within this 
band ll. This band is a measure of the systematic, or intrinsic, sum-rule uncertainty. 



The method of sum rules in QCD contains a set of prescriptions (see e.g. [3]) which are believed to provide such 
a systematic error. In QCD this, however, always remains a conjecture — it is impossible to prove that and even 
to check unambiguously whether the range provided by the standard sum-rule procedures indeed contains the actual 
value of the bound-state parameter. 

The only possibility to check the reliability of the corresponding procedures is to apply the method to a problem 
where the parameters of the ground state may be calculated independently and exactly. Presently, only quantum- 
mechanical potential models provide such a possibility. (For a discussion of many aspects of sum rules in quantum 
mechanics we refer to 0, i, H, 0, 11 H ) ■ 

A simple harmonic-oscillator (HO) potential model is aperfect tool to achieve this goal: it possesses two essential 
features of QCD — confinement and asymptotic freedom [J] — and has the following advantages: (i) the bound-state 
parameters (masses, wave functions, form factors) are known precisely; (ii) direct analogues of the QCD correlators 
may be calculated exactly. 

Making use of this model, we have already studied the extraction of ground-state parameters from different types 
of correlators: namely, the ground-state decay constant from a two-point vacuum-to- vacuum correlator [ll[, the form 
factor from a three-point vacuum-to- vacuum correlator [l2| , and the form factor from a vacuum-to-hadron correlator 
[lU . We have demonstrated that the standard procedures adopted in the literature for obtaining the systematic errors 
do not work properly: for all types of correlators the true known value of the bound-state parameter was shown to 
lie outside the band obtained according to the standard criteria applied. 

One of our most important results was the demonstration of the fact that the "Borel stability criterion" (based 
on a self-evident statement that the physical observables should not depend on the Borel parameter — an auxiliary 
parameter of the method) is a necessary but by far not sufficient criterion to guarantee a good extraction of the 
hadron parameters. Moreover, combined with the assumption of a Borel-parameter independent effective continuum 
threshold it leads to the extraction of very inaccurate values of the bound-state parameters. 

The results reported in [ll], [l2|, [H, [H| gave us a solid ground to claim that also in QCD the actual accuracy of the 
method turns out to be much worse than the accuracy expected on the basis of the standard criteria. 

Notice, however, that these results contained not only the cautious messages concerning the application of the 
method of sum rules to hadron properties, but also the identification of the main origin of the difficulties of the 
method: these difficulties arise from an over-simplified Ansatz for the hadron continuum which is modeled as a 
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perturbative contribution above a Borel-parameter independent effective continuum threshold. We introduced the 
notion of the exact effective continuum threshold, which corresponds to the true bound-state parameters: in a HO 
model the true hadron parameters — decay constant and form factor — are known and thus the exact effective 
continuum thresholds for different correlators may be calculated. We have demonstrated that the exact effective 
continuum threshold (i) is not a universal quantity and depends on the correlator considered (i.e., it is in general 
different for two-point and three-point vacuum-to- vacuum correlators, and for vacuum-to-hadron correlators), and (ii) 
depends on the Borel parameter and, in the case of the form factor, also on the momentum transfer. 

The understanding of this fact allowed us to make a step forward and to propose a new algorithm for extracting 
gro und-state parameters from correlators in those cases where the ground-state energy is known 1 : a simple idea of 
15] is to relax the standard assumption of a Borel-parameter independent effective continuum threshold, and allow 
for a Borel-parameter dependent quantity. The parameters of the T-dependent effective continuum threshold then 
may be found by fitting the average energy of the dual correlator (i.e. the correlator containing only the contribution 
of the low-energy region below the effective continuum threshold) to the known energy of the bound state. In [lj| [l6| 
we have shown that the application of this idea for vacuum-to-vacuum correlators leads to a considerable increase of 
the actual accuracy of the method. 

In the present paper we extend this idea to the extraction of hadron form factors from vacuum-to-hadron correlators 
These correlators are the basic objects for the method of light-cone sum rules, which has been extensively applied 
in the recent years to the analysis of B-decay form factors [ID, [T^] • The issue of the accuracy of light-cone sum rules in 
QCD is extremely important, since the results from this method are extensively used in precision electroweak physics 

We demonstrate that the application of our algorithm yields two crucial improvements: first, it increases significantly 
the accuracy of the extracted form factors; second, it enlarges considerably the range of the momentum transfers where 
the form factor may be extracted from vacuum-to-hadron correlators. 

2. THE MODEL 

Let us recall the essential formulas: We consider a non-relativistic HO model defined by the Hamiltonian (r = \r |) 

H = H + V{r), H =p 2 /2m, V{r) = mw 2 r 2 /2, (1) 
where all characteristics of the bound states are calculable. For instance, for the ground (g) state one finds 

E g = -w, ^ g (r) = {—) exp(-mw 2 r 2 /2), F s (q) = exp(-q 2 /Amu), (2) 

where the elastic form factor of the ground state is defined according to (q = \q\) 

F e (q) = = [ cPklpUk) i> s (k-q), (3) 



with the current operator J(q) given by the kernel 

(f'| J(q)\r) = exp{iq- f) 8 {3) (r-r'). (4) 
Notice that the non-relativistic elastic current is conserved, which leads to the relation F g (q = 0) = 1. 

3. THE VACUUM-TO-HADRON CORRELATOR 

In order to apply the method of sum rules to hadron form factors at intermediate and large monentum transfers a 
Borelized vacuum-to-hadron amplitude of the T-product of two quark currents is used. The analogue of this quantity 
in quantum mechanics has the form [13j 

A(T,q) = (f=0\G(T)J(q)\^). (5) 



1 In our opinion the knowledge of the ground-state mass is mandatory for a successful application of the method of sum rules to bound 
states. Therefore any attempt to study the existence of the ground state with the method of sum rules does not yield trustable results. 



3 



where G(T) = exp(-HT) is the evolution operator in imaginary time T (i.e., the Borel transform E — > T of the full 
Green function of the model, (H — E)~ r ) and J is the current operator defined by ([!]). 

An obvious disadvantage of this correlator compared with vacuum-to-vacuum correlators is the necessity to know 
for its calculation the ground-state wave function. As a bonus, this correlator receives an enhanced contribution of 
the ground state, which makes it potentially more attractive for the extraction of the ground-state form factor. 

The spectral representation for A(T, q) may be written in the form 

A(T,q) = J d 3 kG(k 2 ,T)i/j g (k-q)= [ dza(z,T,q), z = k 2 /2m, (6) 

where H 



G(P,T) = (f= 0\G(T)\k) = ( J )3/2 [cosh( ^ )]3/2 e*P (-^tanh(.T)j . (7) 
Explicit expressions for a(z, T, q) and A(T, q) may then be easily obtained using ip g of Eq. Ip[), For A(T, q) one finds 

The function A(T, q) depends on two dimensionless variables q 2 = q 2 /muj and u>T. 

Expanding in Eq. ([6]) the Green function G(k 2 ,T) in powers of u) generates the analogue of the QCD twist expansion 
for the amplitude A(T,q) (see [l3| for details). 

The ground-state contribution to the correlator has the form 

A s (T,q) = Mr = 0)exp(-E g T)F g (q), (9) 

so one finds 

^'"-^r-^e— V (1 0, 



A(T, q) r V 4tow 

Notice the following features of the correlator A(T, q): 

(i) At large T, the ground state provides the dominant contribution to A, similar to the case of vacuum-to-vacuum 
correlators: 

A(T,q)^A g (T,q), for T -» oo. (11) 

(ii) At small q, because of current conservation, the ground state dominates the correlator for all T. This is a specific 
feature of A which arises due to the choice of the initial state. Thus, compared with vacuum-to-vacuum correlators, 
the correlator A "maximizes" the ground-state contribution. 



4. EXTRACTION OF THE GROUND-STATE FORM FACTOR 



Employing the quark-hadron duality hypothesis, which assumes that the excited-state contribution is dual to the 
high-energy region of quark diagrams, one constructs the correlator dual to the ground state and gets the following 
sum rule for the form factor F g (q): 2 

A dual (T,q,z cS )= ( M{ a(z,T,q)dz = ^ s {r = 0)exp(-E s T)F s (q). (12) 
Jo 

Relation (|12p constitutes the definition of the exact effective continuum thresholds z e s(T, q). 



2 This is the standard but not the unique way to define the dual correlator: another possibility — to apply the cut only to the leading- 
twist contribution and thus to attribute the contributions of all higher twists exclusively to the ground state — has been discussed in 
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Fig. 1: (a) Exact effective continuum threshold z c g(T, q) obtained by solving numerically Eq. (|12|l using the exact bound-state 
parameters ip s (r = 0) and F s (q) vs. Euclidean time T for several values of q = q/^/muj. (b) Relative contribution of the ground 
state to the correlator A s (T,q)/A(T,q) vs. T for several values of q. The vertical lines indicate the window 0.2 < uiT < 0.8. 
Red: 9=1, green: 9 = 2, blue: 9 = 3, light-blue: 9 = 4, black: 9 = 5. 



Let us emphasize the following point: As is clear from Eq. ([6|), the T-dependence of the spectral density a(z, T, q) 
is due to the T-dependence of the evolution operator G(z,T) of Eq. ([7]); the latter may be expanded in a series 
of exponential functions with increasing powers. So, the integration over z of a(z,T,q) in Eq. (|12[) can produce a 
single exponential corresponding to the ground state only if the effective continuum threshold z e g has a nontrivial 
T-dependence. 

More generally and independently of the specific form of the interaction potential, the effective continuum thresh- 
old for the vacuum-to-hadron correlator should depend both on T and q as a consequence of the structure of the 
correlator. This may be seen in two ways: first, the dominance of the ground state in the correlator at small q implies 
z e fi (T, q — > 0) — > oo for all T. Second, for large T, the analytic properties of the dual correlator (which is a hand-made 
nonperturbative object whose properties are rather different from properties of correlators in perturbation theory 
[Til) require that z e g(T — > oo, q) — > oo for all q. 

The explicit T- and (/-dependences of z e g (T, q) can be obtained by solving Eq. (|12[) making use of the exact bound- 
state parameters ip g (r = 0) and F g (q). In the HO model, this can be easily done numerically. The corresponding 
results are shown in Fig. [IJa). Obviously, the effective continuum threshold z c s(T,q) does depend strongly on both 
T and q. 

Let us consider now a restricted problem where the energy E g of the ground state and its wave function at the 
origin, ip s (r = 0), are known, and try to determine its elastic form factor from the sum rule fj 12[) . 

First, according to [lj we should determine the working interval of T — the Borel "window" where the method 
may be applied to the extraction of the ground-state parameter: 

(i) The upper boundary of the T-window is obtained from the condition that the truncated expansion gives a good 
approximation to the exact correlator. Since in the HO model the correlator is known exactly, the upper boundary 
is T = oo. However, to be close to realistic situations, where only a limited number of higher- twist corrections is 
available, following our results reported in [13j we define the upper boundary of the window as loT < 0.8. 

(ii) the lower boundary of this T-window is determined by the condition that the ground state gives a "sizable" 
contribution to the correlator. 3 As can be seen from Fig. Hlb), if we require the ground-state contribution to exceed, 
say, 50%, the window disappears already for q = q/^/mto > 4. We shall see, however, that our algorithm allows one 
to extract the form factor, although with a worse accuracy, also in the region where the ground-state contribution to 
the correlator is below 30%. Thus, our algorithm opens the possibility to study also the region of large momentum 
transfers. 

For practical purposes we need to set a lower boundary of the working region, so we choose the window to be 
0.2<wT<0.8. 

Second, we must choose a criterion to fix the effective continuum threshold z e ff(T, q). We proceed in the following 



3 Notice, however, that in applications of light-cone sum rules in QCD the magnitude of the ground-state contribution is usually not 
considered. 
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Fig. 2: (a) Ratio Fduai(g) / Fg(q) , extracted from the sum rule fl2jl using different approximations z~g(T, q) in Eq. (fT3|) . vs. q. 
Different lines correspond to different Ansatze for the effective continuum threshold: red — constant (n = 0), blue — quadratic 
(n = 2), black — quartic (n = 4). (b) Ratio Fduai(q)/ F g (q), extracted from the sum rule (|12p using different approximations 
z cS (T m Eq. (|13p depending on n, the power of the polynomial used in the fit. The results for two values of the momentum 
transfer are presented: q = 3 — green dots, q = 5 — violet dots. 



way: We consider a set of T-dependent Ansatze for the effective continuum threshold: 

n 

(The standard procedure adopted in all sum-rule applications in QCD is to assume a T-independent quantity.) Now, 
at each value of q we fix the parameters on the r.h.s of (flS")) as follows: we calculate the dual energy 

£dual(T,<z) = - — log^ dual (T,g,^ ff (T, g )), (14) 

for the T-dependent z e g of Eq. (|13p. Then we evaluate -E<j U ai(T, q) at several values of T = T (i = 1, . . . , N, where N 
can be taken arbitrary large) chosen uniformly in the window. Finally, we minimize the squared difference between 
Tduai and the exact value E g : 

1 N 

X 2 = ;tf£[£duai(Ti, 2 - (15) 

8=1 

The results for the dual form factor obtained after optimizing in this way the parameters of z c g are shown in Fig. [21 



A. T-independent effective continuum threshold 

Let us consider first the case of the standard T-independent approximation. One obtains the form factor in the 
region q < 1.5 with better than 10% accuracy, whereas the accuracy falls down rather fast with increasing q (see 
Fig. (Ha)). The real problem is, however, that the magnitude of the error cannot be guessed on the basis of the 
standard criteria adopted in the method: e.g., at q = 2, the variation of the extracted form factor in the window 
is only about 2%, which mimics an accurate extraction of the form factor, whereas the actual error comprises 15%. 
Thus we conclude that, similar to the case of vacuum-to-vacuum correlators, the error of the form factor extracted 
from the vacuum-to-hadron correlator cannot be determined from the variation of the dual form factor in the window. 



B. T-dependent effective continuum threshold 

We are going to discuss now the quadratic and quartic Ansatze, zj%! (T, q) and z^g (T, q) (the results obtained for 
n = 1 are close to the results for n — 0, and the results for n — 3 are close to those for n = 2). As can be seen from 
Fig. (Ha), the form factor may be extracted with reasonable accuracy (better than 15%) up to q — 5. This relatively 
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high accuracy is reached in spite of the facts that (i) the relative contribution of the ground state to the correlator 
falls down considerably below 50% in the T-window, and (ii) the form factor at q = 5 falls down compared to its value 
at q = by almost three orders of magnitude. 

It should be mentioned that both for n = 2 and n = 4 the dual energy E^ ua i(T,q) and the dual form factor 
Fdu&\(T,q) are extremely stable (better than 0.1%) in the window. In spite of this stability, the actual error of the 
extracted form factor for both n = 2 and n = 4 is found to be at the level of 10-15%. Again, this error could not 
be guessed on the basis of the standard Borel stability criterion, which thus does not work for vacuum-to-hadron 
correlators. Previously, we made a similar observation for the case of vacuum-to- vacuum correlators (Til ll~2l [T5| . 
Interestingly, the actual error may be probed if one compares the results obtained for quadratic and quartic Ansatze 
for z c g] moreover, in the case under consideration these results provide the actual band of values which contains the 
true hadron form factor. 

Fig. HJb) shows the results for the form factor obtained by the fit as function of the power n of the T-polynomial 
in ([13)) for two values q = 3 and q — 5. Let us notice that the fitting procedure becomes unstable for n > 7: the fit 

leads to coefficients Zj in (I13|) which start to increase with j. Respectively, increasing the power of the polynomial 
does not lead to an improvement of the accuracy of the extracted from factor. This behaviour might be expected and 
reflects the fact that the extraction of the form factor from the correlator known in the limited range of T cannot be 
done arbitrarily accurate. 

5. CONCLUSIONS 

We discussed the extraction of the form factor from the vacuum-to-bound-state correlator, making use of a quantum- 
mechanical potential model. We formulated a new algorithm for fixing the effective continuum threshold — one of 
the basic ingredients of the method of sum rules — which to a great extent determines the numerical value of the 
form factor obtained by this method. 

The main messages of our analysis are as follows: 

• We would like to emphasize that the exact effective continuum threshold in vacuum-to-bound-state correlators, 
(as well as the exact continuum threshold in vacuum-to-vacuum correlators) depends on the Borel parameter T. This 
conclusion is based on both the general properties of the vacuum-to-hadron correlator and on explicit calculations. 

• Assuming a T- independent (i.e., a Borel-parameter independent) Ansatz for the effective continuum threshold, 
the error of the ground-state form factor extracted from the vacuum-to-hadron correlator turns out to be typically 
much larger than (i) the error of the description of the exact correlator by the truncated twist expansion and (ii) the 
variation of the dual form factor in the Borel window. Combining this statement with our previous results obtained 
for vacuum-to- vacuum correlators [TTl [T^, [H, [l6j], we conclude that in all versions of sum rules "Borel stability" 
cannot be used as a criterion of the reliability of the extracted form factor. 

• Allowing for a T-dependent effective continuum threshold and fixing it according to Eq. ([15")) leads to dramatic 
improvements in the extracted ground-state form factor: 

First, it extends considerably the range of the momentum transfers where the form factor may be extracted with 
an accuracy at the level of 10-15%. In particular, it extends the applicability of the method also to the region of large 
momentum transfers, where the ground state does not give the dominant contribution to the correlator. 

Second, it improves considerably the accuracy and reliability of the form-factor extraction compared to the standard 
T-independent Ansatz for the effective continuum threshold that has been used in all previous applications of sum 
rules in QCD. Moreover, in the example under consideration, we managed to obtain the band of values which contains 
the true form factor, by analysing different Ansatze for the effective continuum threshold. If this property holds in 
QCD, it would eventually provide a means to gain a rigorous control over the accuracy of light-cone sum rules. This 
important issue requires further investigation. 

The application of the proposed ideas to vacuum-to-hadron correlators in QCD seems very promising. We believe 
that it will lead to a considerable improvement of hadron form factors obtained from light-cone sum rules. 
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